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1.  OBJECTIVE 


The  goal  of  this  study  is  to  develop  accurate  computational  techniques  for 
modeling  acoustic  scattering  from  a  target  embedded  in  a  two-phase  poroelastic 
medium,  such  as  a  water-filled  sandy  sediment.  Such  modeling  capabilities 
would  substantially  enhance  existing  sonar  performance  prediction  models. 

Previous  tasks,  ending  with  TD  No.  01A1038,  demonstrated  the 
penetration  of  significant  acoustic  energy  into  sandy  sediments  at  shallow 
grazing  angles.  This  phenomenon  has  potentially  important  consequences  when 
applied  to  the  detection  and  classification  of  buried  targets.  It  opens  the  door  to 
the  possibility  of  using  active  sonar  backscattering  to  detect  and  categorize 
buried  objects  at  long  range.  However,  this  possibility  can  only  be  realized  if  the 
capacity  to  accurately  model  acoustic  penetration  and  scattering  is  well 
developed. 

Contrary  to  experimental  observations,  current  elastic  medium 
propagation  models  such  as  SAFARI  and  PE  are  unable  to  account  for  acoustic 
penetration  at  shallow  grazing  angles.  This  is  an  inherent  limitation  of  elastic 
sediment  models.  For  a  given  frequency,  these  models  predict  the  existence  of 
one,  and  only  one,  acoustic  compression  wave  within  the  sediment.  For  smooth 
surfaces,  model  parameters  that  are  compatible  with  observed  values  of 
sediment  acoustic  wave  speed  and  attenuation  inevitably  predict  a  critical  grazing 
angle,  below  which  negligible  acoustic  energy  penetrates  the  sediment. 

Experimental  measurements  also  indicate  the  existence  of  a  second 
acoustic  wave  in  water-saturated  sand.  Therefore,  a  model  that  can 
accommodate  two  acoustic  waves  is  required.  One  such  model  is  Biot's  theory 
of  acoustic  propagation  in  poroelastic  media,  which  predicts  the  existence  of  both 
a  fast  and  a  slow  acoustic  compression  wave.  Model  parameters  that  match  the 
fast  wave  speed  with  observed  values  predict  a  critical  grazing  angle  for  the  fast 
wave,  but  not  necessarily  for  the  slow  wave.  Thus,  acoustic  energy  can 
penetrate  into  the  sediment  at  any  grazing  angle. 

This  study  explores  the  effectiveness  of  Biot's  theory  in  modeling  acoustic 
scattering  from  spherical  targets  embedded  in  a  sandy  sediment.  If  extended 
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testing  proves  Biot’s  theory  to  be  reliable,  the  ultimate  goal  is  to  combine  this 
scattering  model  with  an  acoustic  propagation  and  penetration  model,  producing 
a  fully  integrated  target  detection  scheme  for  targets  buried  in  sandy  sediments. 
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2.  INTRODUCTION 


Over  the  years,  the  ocean  sediment  has  been  modeled  as  a  fluid,  elastic 
solid,  and  most  recently  as  a  fluid-filled  poroelastic  solid. ''•3  Although  these 
models  make  similar  predictions  in  some  respects,  there  are  notable  differences. 
Perhaps  the  most  significant  is  the  number  and  types  of  wave  disturbances  each 
type  of  medium  will  support.  Fluid  media  support  a  single  compression  wave, 
elastic  media  support  one  compression  wave  and  one  shear  wave,  and  poro¬ 
elastic  media  support  two  compression  waves  and  one  shear  wave.  The  indices 
of  refraction  and  scattering  properties  of  these  different  waves  can  vary 
considerably.  The  three  models  result  in  decidedly  different  scenarios  for  the 
detection  of  objects  buried  within  the  seabed.  This  is  especially  true  for  waves 
incident  on  a  smooth  seafloor  at  shallow  grazing  angles. 

Recent  analysis  of  experimental  bottom  penetration  results  by  Chotiros^’® 
is  consistent  with  the  poroelastic  theory  developed  by  Biot.^  ^  j^us  jt  seems 
reasonable  to  apply  Biot's  theory  to  model  scattering  from  an  inclusion  within  the 
sediment.  Indeed,  a  few  approaches  to  this  problem  have  already  been 
investigated.  In  the  long  wavelength  limit,  analytic  expressions  were  obtained  by 
Berryman^  for  scattering  from  a  spherical  elastic  inhomogeneity  within  a  Biot 
medium.  More  recently,  KargI  and  Limio  developed  a  transition  matrix  technique 
and  applied  it  numerically  to  scattering  from  a  spherical  Biot  inhomogeneity  within 
a  Biot  medium.  But  the  most  general  treatment  to  date  was  given  by  Zimmerman 
and  Stern,'' I  who  derived  a  series  solution  for  the  scattering  of  plane 
compression  waves  by  a  fluid,  elastic,  or  rigid  spherical  inclusion  within  a  Biot 
medium. 

The  numerical  results  presented  here  are  based  on  the  work  of 
Zimmerman  and  Stern.  The  incident  wave  can  be  composed  of  any  linear 
combination  of  Biot  fast  and  slow  compression  waves.  Also,  the  scattered  wave 
is  easily  decomposed  into  its  fast,  slow,  and  shear  components.  The  angular 
dependence  of  the  scattered  fast  and  slow  components  is  investigated,  and 
compared  with  experimental  results. 
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3.  MODEL 


Biot's  theory  predicts  acoustic  propagation  through  a  poroelastic  medium 
as  a  coupled  wave  motion  within  the  solid  and  the  pore  fluid.  The  medium  is 
modeled  as  a  solid  with  tubular  pores,  as  illustrated  in  Fig.  3.1.  The  propagating 
wave,  in  the  direction  of  the  pore  tubes,  can  be  decomposed  into  three 
components  --  fast  and  slow  compression  waves,  and  a  shear  wave.  The 
equations  of  motion,  in  a  form  used  by  Zimmerman  and  Stern, '■  are 

v.<T, +  B,  V)  =  ^(p;,u  +  p,,Vl  (3.1) 

OC  at 

and 

+  B.  +  Z? -^(U  -  V)  =  -^(p;,u  +  p;.  VI.  (3.2) 

di '  dr 

where  U  and  V  are  the  solid  and  fluid  displacements,  respectively,  Bg  and  B/^are 
the  body  forces  per  unit  volume  of  bulk  material  applied  to  the  solid  and  fluid, 
respectively,  b  \s  a  viscous  coupling  parameter  relating  drag  force  to  relative 
solid-fluid  velocity,  s  is  the  negative  pore  pressure  times  the  porosity,  Os  is  the 
solid  stress  tensor,  and  the  pS  are  mass  coefficients. 

This  work  is  a  numerical  study  of  the  scattering  of  transient  plane 
compression  waves  from  an  elastic  sphere  embedded  in  an  unbounded  Biot 
medium.  The  geometry  of  the  orobiem  is  shown  in  Fig.  3.2.  The  incident  wave 
can  be  any  linear  combination  of  Biot  fast  and  slow  waves.  A  fast  Fourier 
transform  is  applied  to  the  incident  wave,  and  each  frequency  component  is 
transformed  to  spherical  coorainates  via  a  Legendre  polynomial  /  spherical 
Bessel  function  decomposition. '‘2  Then  the  steady-state  solution  of  the  above 
equations  of  motion  are  determined  for  each  frequency  by  applying  the 
appropriate  boundary  conditions  at  the  Biot-sphere  interface.  The  details  of  this 
solution  technique  were  derived  by  Zimmerman  and  Stern,'' ^  using  an  approach 
previously  applied  by  Ying  ana  Trueihs  and  Pao  and  Mow''*  *or  the  case  of  a 
spherical  object  embedded  in  an  elastic  medium.  The  steady  state 
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wave  wave 
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Figure  3.1 

Illustration  of  Biot's  theory  of  acoustic  propagation 
through  a  poroelastic  medium. 
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Spherical  Inclusion 


Figure  3.2 
Problem  geometry. 


ARL  -  UT 
AS-95-1235 
DJY/RFO 
12-95/svr 


solution  for  each  angular  frequency,  con,  is  multiplied  by  exp(-/cjnt),  and  the  results 
are  summed  to  give  the  full  transient  scattered  wave  solution. 
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4.  NUMERICAL  IMPLEMENTATION 


The  numerical  implementation  of  the  series  solution  developed  by 
Zimmerman  and  Stern  was  quite  straightforward,  although  somewhat  tedious. 
The  only  serious  difficulties  that  could  potentially  compromise  the  calculations 
were  the  choice  of  Bessel  function  and  Legendre  polynomial  routines,  and  the 
choice  of  a  convergence  criterion  for  the  series  solution.  The  reliability  of  the 
numerical  results  critically  depends  on  the  reliability  of  these  two  choices. 

The  Bessel  function  and  Legendre  polynomial  routines  used  were  the 
library  of  routines  developed  by  Amos,  of  Sandia  National  Laboratories.  The 
routines  are  in  double  precision  and  were  used  to  evaluate  Bessel  functions  of 
the  first,  second,  and  third  kind,  as  well  as  Legendre  polynomials. 
The  routines  are  available  on  the  World  Wide  Web  at  the  URL 
"http://netlib.att. com/magic/netlib_find?db=0&pat=amos.‘'  A  general  overview  of 
the  routines  and  a  discussion  of  their  capabilities  and  limitations  are  given  in  the 
file  "amos/readme." 

The  solution  series  was  considered  to  have  converged  when  both  the 
incremental  radial  solid  displacement  term  and  the  incremental  radial  stress  term 
had  decreased  in  magnitude  to  a  point  less  that  10'^  times  the  series  sum  prior  to 
those  terms.  For  a  3.5  cm  diam  sphere,  and  in  the  frequency  range  from 
100  kHz  to  300  kHz,  convergence  was  obtained  out  to  a  distance  of  four  sphere 
diameters  by  using  30  terms  for  the  incident  wave  series  and  ten  terms  for  the 
scattered  wave  series.  All  the  computations  were  done  in  double  precision. 
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5.  NUMERICAL  RESULTS 


The  scattering  of  transient  plane  compression  waves  by  a  3.5  cm  diam 
steel  sphere  embedded  in  an  unbounded  Biot  medium  was  simulated 
numerically.  The  properties  of  the  Biot  medium  and  the  steel  sphere  are  given  in 
Tables  5.1  and  5.2,  respectively.  The  input  was  a  vibrating  planar  surface  at  a 
distance  of  10  cm  from  the  sphere.  The  Fourier  transform  of  the  input  vibration 
was  modeled  as  a  cosine  half-wave  in  the  frequency  domain,  in  the  frequency 
band  100-300  kHz,  as  shown  in  Fig.  5.1.  Of  the  incident  plane  wave,  98%  of  the 
energy  incident  on  the  sphere  was  in  the  fast  wave  component,  2%  in  the  slow 
wave  component,  and  a  negligible  amount  in  the  shear  wave  component.  The 
angular  dependence,  as  a  function  of  time,  of  the  scattered  fast  and  slow 
compression  waves  was  investigated.  The  results  are  shown  in  Figs.  5.2 
and  5.3. 

For  the  scattered  fast  wave,  the  specular  return  and  the  successive 
Rayleigh  waves  were  found  to  decrease  as  the  angle  changed  from  180° 
(backscatter)  to  90°.  As  the  angle  changed  from  90°  to  0°  (forward  scatter),  the 
specular  return  increased  again,  whereas  the  behavior  of  the  Rayleigh  waves 
was  unclear.  For  all  angles,  most  of  the  scattered  energy  was  found  to  be  in  the 
specular  return.  The  time  difference  between  the  arrival  times  of  the  specular 
return  at  different  angles  was  consistent  with  refraction  of  the  plane  wave  through 
the  sphere  with  a  longitudinal  wave  speed  of  6001  m/s.  Likewise,  the  time 
difference  between  the  specular  return  and  the  successive  Rayleigh  waves  was 
consistent  with  a  transverse  wave  moving  around  the  circumference  of  the 
sphere  with  a  wave  speed  of  3208  m/s. 

For  the  scattered  slow  wave,  the  specular  return  was  found  to  increase  as 
the  angle  changed  from  180°  to  90°,  whereas  the  Rayleigh  waves  were  found  to 
decrease.  From  180°  to  135°,  most  of  the  scattered  energy  was  in  the  first 
Rayleigh  wave,  but  at  all  other  angles  it  was  mostly  in  the  specular  return.  As  the 
angle  changed  from  90°  to  0°,  the  specular  return  decreased  and  the  first 
Rayleigh  wave  decreased.  The  behavior  of  the  other  Rayleigh  waves  in  this 
regime  was  unclear.  The  time  difference  between  the  arrival  time  of  the  specular 
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Table  5.1 

Material  properties  of  the  simulated  Biot  medium. 


Fluid  viscosity 

1.00x10-3  kg/ms 

Fluid  mass  density 

1 000  kg/m3 

Fluid  bulk  modulus 

2.25x109  Pa 

Grain  mass  density 

2650  kg/m3 

Grain  bulk  modulus 

7.00x109  Pa 

Frame  shear  modulus 

2.61  X  107  Pa 

Frame  bulk  modulus 

5.30  X  109  Pa 

Frame  shear  attenuation 

1 .303  dB/wavelength 

Frame  bulk  attenuation 

1 .303  dB/wavelength 

Frame  porosity 

0.44 

Frame  virtual  mass  parameter 

1.636 

Frame  pore  size  parameter 

3.209  X  10-4  m 

Frame  permeability 

2.265  X  10-9  m2 

Table  5.2 

Properties  of  the  simulated  steel  sphere. 


Diameter 

3.5  X  10-2  m 

Mass  density 

7850  kg/m3 

Longitudinal  wave  speed 

Transverse  wave  speed 

6001  m/s 

3208  m/s 
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(a)  Time  aomain 


(b)  Frequency  aomain 


Figure  5.1 

Incident  wave  at  Zq  =  -10  cm. 
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return  at  different  angles  was  consistent  with  refraction  of  the  plane  wave  through 
the  sphere  with  a  longitudinal  wave  speed  of  6001  m/s.  Likewise,  the  time 
difference  between  the  specular  return  and  the  successive  Rayleigh  waves  was 
consistent  with  a  wave  moving  around  the  circumference  of  the  sphere  with  a 
transverse  wave  speed  of  3208  m/s. 
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6.  EXPERIMENTAL  RESULTS 


At  present,  experimental  comparisons  for  the  full  angular  spread  of 
predictions  given  in  Figs.  5.1  and  5.2  are  not  available.  However,  both 
experimental  measurements  and  numerical  data  were  obtained  for  the  forward 
scattering  case  (0°). 

The  experiments  were  conducted  in  a  glass  aquarium  filled  with  sand  and 
water.  The  aquarium  was  25  cm  wide,  51  cm  long,  and  30  cm  deep.  It  was  filled 
to  a  depth  of  20  cm  with  a  sand  and  water  mixture,  with  an  additional  6  cm  of 
water  above.  Trapped  air  bubbles  were  coaxed  from  the  sand  using  water 
circulation  pumps  and  a  vacuum  chamber.  With  the  circulation  pumps  in  the 
water,  the  aquarium  was  placed  in  the  vacuum  chamber  and  the  pressure  was 
reduced  to  0.06  atm.  The  circulation  pumps  were  then  turned  on  for  3-4  h, 
enhancing  the  coalescence  and  removal  of  air  bubbles.  Afterwards,  the  pumps 
were  turned  off  and  the  aquarium  remained  under  vacuum  overnight,  allowing 
more  trapped  air  to  escape.  Next,  the  aquarium  was  removed  from  the  vacuum 
chamber,  and  the  source,  target,  and  receiver  assembly  were  inserted  into  the 
sand  very  carefully,  in  the  configuration  shown  in  Fig.  6.1.  The  sand  was  then 
encouraged  to  settle  by  placing  the  filled  aquarium  on  a  platform  vibrating  with  a 
very  small  vertical  amplitude  at  a  frequency  of  58  Hz  for  10  min.  The  aquarium 
and  vibrating  platform  are  shown  in  Fig.  6.2.  Finally,  the  entire  apparatus  was 
placed  back  in  the  vacuum  chamber,  and  the  pressure  was  reduced  to  0.06  atm. 
Measurements  were  not  taken  until  the  apparatus  had  been  under  vacuum  for  an 
additional  2  h,  allowing  more  air  bubbles  to  escape  from  the  sediment. 

The  total  assembly,  shown  in  Fig.  6.3,  was  designed  to  maintain  equal 
depths  for  each  of  the  three  elements.  The  steel  rod  supporting  the  target  was 
0.6  cm  in  diameter  at  the  rod/target  interface.  The  steel  sphere  target  was 
3.5  cm  in  diameter.  The  sound  source  and  receiver  are  wideband  transducers 
constructed  of  a  new  composite  material  containing  piezoelectric  particles  in  a 
rubber  matrix.  The  assembly  was  engineered  so  that  the  distances  between  the 
source  and  target  and  between  the  target  and  receiver  were  adjustable,  as  was 
the  angle  between  the  source,  target,  and  receiver.  It  is  capable  of  measuring 
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Figure  6.1 

Photograph  of  experimental  source,  target,  and  receiver  assembly. 
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Steel  supporting  rods 


Figure  6.3 

A  cross-sectional  schematic  of  the  experimental  apparatus. 
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scattered  signals  as  a  function  of  the  bistatic  angle  d,  as  shown  in  Fig.  6.4.  The 
experiment  reported  here  used  9  =  0°. 

Two  sets  of  measurements  were  taken  --  one  with  a  spherical  target,  in 
water-saturated  sand,  and  the  other  with  no  targer.  in  water.  In  both  cases,  the 
distance  between  the  faces  of  the  source  and  receiver  transducers  was  20  cm. 
When  used,  the  target  was  midway  between  the  two  transducers.  The  measured 
data  consisted  of  the  input  current  of  the  source  transducer  and  the  output 
voltage  of  the  receiver  transducer.  The  data  were  processed  by  first  filtering  out 
the  (very  small)  dc  component  and  then  tapering  the  trailing  end  of  the  signal  to 
compensate  for  periodic  mismatch.  These  processed  data  were  then  Fourier 
transformed,  and  the  transfer  functions  were  calculated  via  a  Wiener  filter. 

The  transfer  function  for  the  case  of  no  target  in  water  was  used  to 
deconvolve  the  transducer  response  from  the  oata  for  the  case  of  a  spherical 
target  in  water-saturated  sand.  The  process  useo  assumes  the  system  is  linear. 
The  transfer  function  for  the  no-target  case  represents  the  unobstructea 
response  of  the  transducers  in  an  essentially  lossless  medium  (water),  relative  to 
water-saturated  sand.  However,  it  does  incluoe  a  spreading  loss  which  is 
inversely  proportional  to  the  square  of  the  separation  between  the  two 
transducers.  This  loss  was  divided  out  of  the  no-iarget  transfer  function,  giving 
the  transfer  function  due  solely  to  the  transducer  response.  This,  in  turn,  was 
divided  out  of  the  transfer  function  for  the  case  of  a  spherical  target  in  water- 
saturated  sand,  via  a  Wiener  filter.  The  result,  assuming  the  system  is  linear, 
was  the  experimental  transfer  function  from  input  cressure  to  output  pressure  for 
a  spherical  steel  target  in  water  saturated  sand. 

The  same  transfer  function  was  simulatea  numerically,  and  the  power 
spectrum  of  the  simulation  was  compared  to  the  experimental  power  spectrum. 
To  account  for  the  narrower  frequency  range  of  the  simulated  incident  wave,  born 
transfer  functions  v/ere  passed  through  a  50-250  kHz  raised  cosine  filter  prior  to 
calculating  the  power  spectra.  The  result  is  shown  in  Fig.  6.5. 

As  shown  by  the  difference  spectrum  in  Fig.  6.5,  the  two  power  spectra 
agree  to  within  ±10  dB  at  all  frequencies.  However,  the  structure  of  the  two 
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Figure  6.4 

Plan  view  of  target,  source,  and  receiver  geometry. 
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Figure  6.5 

Comparison  of  the  experimental  and  simulated  power  spectral 
responses  in  the  case  of  forward  scattering  (0°). 
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spectra  are  notably  different.  The  simulated  spectrum  has  prominent,  closely 
spaced  peaks  at  83  kHz  and  107  kHz,  and  a  third  prominent  peak  at  207  kHz, 
whereas  the  experimental  spectrum  only  has  two  large  bumps.  The  bumps  do, 
however,  correspond  roughly  to  the  peaks  in  the  simulated  spectrum.  The  larger 
bump  overlaps  the  two  closely  spaced  peaks,  and  the  smaller  bump  overlaps  the 
third  peak. 

A  comparison  of  the  bandlimited  impulse  responses  of.theoretical  model 
and  experiment  is  shown  in  Fig.  6.6.  It  indicates  that  there  is  much  similarity  in 
the  initial  part  of  scattereo  pulse.  The  main  difference  is  that  the  model  predicts  a 
more  persistent  ringing  than  the  experimental  data,  indicating  that  the 
experimental  system  is  more  heavily  damped.  Possible  reasons  include 
(1)  imperfections  in  the  coupling  between  tangential  motion  at  the  surface  of  the 
sphere  and  the  surrounding  sand,  resulting  in  greater  damping  in  the  Rayleigh 
wave  response  than  preoicted,  and  (2)  distortion  causeo  by  the  weld  between 
sphere  and  mounting  rod. 
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Time  -  us 


Figure  6.6 

Comparison  of  the  experimental  and  simulated  impulse 
responses  in  the  case  of  forward  scattering  (0°). 
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7.  CONCLUSIONS 


The  results  in  Figs.  5.2  and  5.3  demonstrate  very  significant  differences  in 
the  behavior  of  the  scattered  Biot  fast  and  slow  waves.  This  diversity  guarantees 
there  will  be  notable  discrepancies  between  the  scattering  predictions  of  a  Biot 
sediment  model  and  an  elastic  sediment  model,  since  the  latter  is  inherently 
incapable  of  manifesting  such  diversity.  If  borne  out  experimentally,  the  strongly 
backscattered  first  Rayleigh  wave  in  the  slow  wave  spectrum  could  be  a 
significant  resource  for  the  long  range  detection  of  targets  buried  in  the  seabed. 
By  the  principle  of  reciprocity,  this  signal  return  would  register  most  strongly  at 
the  source  of  the  incident  wave.  I.e.,  at  the  location  of  the  active  sonar  used  to 
scan  the  sediment.  This  is  the  most  advantageous  condition  for  self-contained 
target  detection.  Further  simulations  and  experiments  are  required  to  determine 
the  relative  differential  scattering  cross  sections  for  incident  fast  versus  slow 
waves. 

The  forward  scattering  simulation  result,  depicted  in  Figs.  6.5  and  6.6,  are 
in  good  qualitative  agreemeni  with  the  experimental  measurements.  The 
simulated  signal  strength  is  closely  centered  on  the  experimental  signal  strength, 
and  the  coarse  grain  features  of  the  two  spectra  are  similar.  However,  there  are 
notable  differences  on  a  finer  scale.  The  two  closely  spaced  peaks  in  the 
simulated  power  spectrum  were  not  resolved  experimentally.  The  reason  for  this 
is  unclear,  but  there  are  several  possibilities.  One  is  imperfections  in  the 
coupling  between  tangential  motion  at  the  surface  of  the  sphere  and  the 
surrounding  sand,  resulting  in  greater  damping  in  the  Rayleigh  wave  response 
than  predicted.  Another  is  that  the  target  support  rod  interfered  with  the 
experimental  resolution  of  the  two  peaks.  Further  experimentation  is  required  to 
determine  whether  one  of  these,  or  some  other,  explanation  is  correct. 
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